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ABSTRACT 

Context. Scientific exploitation of the ever increasing volumes of astronomical data requires efficient and practical methods for data 
access, visualisation, and analysis. Hierarchical sky tessellation techniques enable a multi-resolution approach to organising data on 
angular scales from the full sky down to the individual image pixels. 

Aims. We aim to show that the Hierarchical progressive survey (HiPS) scheme for describing astronomical images, source catalogues, 
and three-dimensional data cubes is a practical solution to managing large volumes of heterogeneous data and that it enables a new 
level of scientific interoperability across large collections of data of these different data types. 

Methods. HiPS uses the HEALPix tessellation of the sphere to define a hierarchical tile and pixel structure to describe and organise 
astronomical data. HiPS is designed to conserve the scientific properties of the data alongside both visualisation considerations and 
emphasis on the ease of implementation. We describe the development of HiPS to manage a large number of diverse image surveys, 
as well as the extension of hierarchical image systems to cube and catalogue data. We demonstrate the interoperability of HiPS and 
Multi-Order Coverage (MOC) maps and highlight the HiPS mechanism to provide links to the original data. 

Results. Hierarchical progressive surveys have been generated by various data centres and groups for ~200 data collections including 
many wide area sky surveys, and archives of pointed observations. These can be accessed and visualised in Aladin, Aladin Lite, and 
other applications. HiPS provides a basis for further innovations in the use of hierarchical data structures to facilitate the description 
and statistical analysis of large astronomical data sets. 

Key words. Surveys - Atlases - Astronomical databases - Catalogs - Virtual observatory tools 


1. Introduction 


Astronomical surveys cover the sky in a rich tapestry of data as the recorded signals of photons with spectral and temporal char¬ 
acteristics, stitched together into multi-dimensional maps for science. Large surveys may cover the entire sky, others span irregular 
patchwork patterns of coverage, and the collections of pointed observations in observatory archives typically project sparse sky 
coverage maps composed of thousands of instrumental fields of view. Systems for managing and manipulating astronomical data 
that cover this wide range of spatial scales from the individual pixels up to the full sky, which may also have spectral or temporal 
axes, are necessary for astronomical research that uses data from many surveys and archives. The sheer multitude of surveys and 
archived data requires that these systems are efficient, interoperable, and easy to implement and use. 

Current and planned astronomical projects are ushering astronomy into an era of petabyte surveys, where source catalogues may 
contain 10 9 -10 12 entries. To highlight a few examples: The Gaia mission ( Perryman et al.|[200l| > is taking data and will measure 
the positions of about one billion stars; LOFAR (van Haarlem et ak] 2013| l"science data products are expected to grow at five 
petabytes yr~ 1 and will generate petabyte size catalogues; the EuclitQmission complete survey will comprise hundreds of thousands 
of images (several tens of petabytes of data), and the expected catalogues will have some ten billion sources. LSST (Ivezic et ah 


2008) will revolutionise the exploration of the transient sky by accumulating some 32xl0 12 photometric measurements, and the 


SKA^is planning surveys with exabyte volumes. Scientific exploitation of these rich surveys, in particular the pursuit of research 
that is based on their combination and comparison, requires that these data are highly interoperable and that the systems that enable 
access to these data be tightly integrated with the whole process of data collection, processing, curation and archiving, and analysis. 


* Corresponding authors e-mail: Mark.Allen@astro.unistra.fr, Pierre.Fernique@astro.unistra.fr 

1 http://sci.esa.int/euclid/ 

2 https://www.skatelescope.org 
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Fig. 1 . HEALPix tessellation of the sphere. Each of the four orthographic projections of the sphere shown here display the twelve base resolution 
HEALPix quadrilaterals shaded with different grey levels. The left most frame shows how the HEALPix tessellation relates to a generic spherical 
longitude and latitude coordinate system. The second, third, and forth spheres to the right are overlaid with HEALPix grids of orders k = 0,1,2 
respectively, to illustrate the hierarchical structure of HEALPix where each pixel is divided into four self-similar pixels at each successive order. 


There are already hundreds of large surveys available though astronomy data centres and services. It is also becoming common 
to consider the observational archives of the pointed observations of a given instrument as a kind of heterogeneous survey data. 
This is also true of the source catalogues derived from cumulated archives of pointed observations, for example the Chandra source 
catalogue ( {Evans et al.|2010 ) and the Hubble source catalogue ( Budavari & Lubow|2012) >. With increasingly interoperable data in 
observational archives, an ultimate goal could involve characterising observational data along many physical axes so that the data 
may be described within a generalised multi-dimensional parameter space, and data could be combined along any compatible axis. 
Presently, there is growing awareness of the potential of practical schemes for characterising the spatial sky coverage of the wide 
diversity of astronomical survey and pointed observation archival data. Our position is that characterising the spatial coverage of 
this diversity of data using a hierarchical scheme leads to innovative solutions to accessing, visualising, and manipulating the data, 
and allows for general cross-comparison of data that scales to many hundreds of data sets. 

Interoperability of astronomical image data is already very good due to compatible formats (FITS) and the use of standardised 
astronomical coordinate systems (Pence et al. |2010| ). Many surveys and archives are made available through web interfaces with 
various kinds of database systems, and virtual observatory (VC^]) metadata and protocols are helping to standardise these interfaces, 
leading to services that provide coordinated access to multiple surveys and data sets. These interfaces typically rely on queries 
expressed in astronomical coordinates and use a query, select, download pattern for accessing cut-out regions of the data. For large 
area studies this approach can be cumbersome because downloading the data for large sky regions at the original pixel resolution 
is inefficient and it is difficult to use a survey as a whole, let alone multiple surveys at the same time. Visualisation interfaces 
have however made significant progress towards interactive displays of multiple surveys across the entire sky. The hierarchical 
techniques that have been employed to enable fast multi-resolution visualisation of astronomy data can be extended to enable easy 
publication and access to image survey data. These techniques can also be generalised to describe heterogeneous archival data, 
multi-dimensional cube data, as well as catalogue data, and have the potential to provide a new level of interoperability between 
these different data types. 

Over the past ten years there has been much improvement in the tools used for the visualisation of astronomical data (Hassan 


& Fluke||2011 >. For large astronomical imaging surveys, it is now commonplace to browse, pan, and zoom into the survey data 


using a tool that accesses the data over the internet and streams only the portions of the data that are needed for the current user 
view. This kind of remote visualisation enables the exploration of large data sets from the full wide scale view of an entire survey 
displayed on a grid of the whole sky, down to the detailed zoomed-in view at the finest spatial limits of the images. Examples of 
such interactive interfaces include Google sky|^] Microsoft World Wide Telescop^j(WWT), Sky-Map/Wiki sk}|^] an d Aladirjj These 
tools and interfaces are made available as client software, in-browser tools, and most recently as applications for mobile devices. 
These interfaces serve a range of visualisation purposes for professional research astronomy, as well as education, citizen science, 
and public outreach, with their ease of use being a major factor in their up-take. 

All of the sky browsing tools mentioned above must balance the considerations of internet bandwidth and the performance of 
the client displays so that the real-time update of the display for panning and zooming is fast enough to be interactive. Important 
considerations in the design of systems for interactive visualisation of remotely stored data are the underlying organisation of the 
data, and how the sky is divided up into regions. All of the applications mentioned above employ a hierarchical multi-resolution 
tessellation of the sky. Several different schemes are in u se: Aladin uses the Hierarchical Equal Area isoLatit ude Pixelization sch eme, 
HEALPit^ (Gorski et al. 2005 |Gorski & Hivon 201 1| ; WWT uses a hierarchical triangular mesh (HTM) (Kunszt et al. 2001 1 and 
an extension of this system called tessellated octahedral adaptive subdivision transform (TOAST); Google sky uses a cylindrical 
method. These different tessellations of the sky have been developed for different purposes and some aspects of their relative merits 
are discussed in O’Mullane et al. (2001). For sky browsing applications, the choice of tessellation scheme dictates how the original 
data must be mapped onto a set of hierarchically organised tiles, and determines the spatial index to be used in the applications. 


3 http://www.ivoa.net 

4 http://www.google.com/sky 

5 http://www.worldwidetelescope.org 

6 http://sky-map.org 

7 http://aladin.unistra.fr 

8 http://healpix.sourceforge.net 
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This paper describes the Hierarchical progressive survey (HiPS) system that has been designed to facilitate the access and 
visualisation of astronomical survey data. HiPS is based on the HEALPix sky tessellation, and is essentially a mapping of survey data 
at various spatial resolutions into a collection of HEALPix tiles. HiPS has been developed at the Centre de donnees astronomiques 
de Strasbourg (CD50I to support the visualisation of all sky imaging survey data in Aladin (Fernique et al. 2010), and has been 
extended to also support source catalogues and three-dimensional cube data. HiPS is designed specifically for astronomical data in 
that it is fundamentally spherical, it takes the astrometric and photometric properties of the original data into account, and emphasis 
is placed on ease of use with no need for special servers or database systems. Section[2]of this paper specifies the hierarchical spatial 
index scheme and its implementation as a file system directory and file structure. Section[3]describes the relationship between HiPS 
and the Multi-Order Coverage (MOC; Fernique et al. 2014[ ) maps that define regions of the sky. Section [4] describes how HiPS is 
applied to astronomical images, catalogues, and three-dimensional data cubes. Section[5]highlights the ~200 HiPS data sets that are 
currently available from various astronomy data centres and outlines how these can be used via the Aladin, Aladin Lite, and other 
compatible client applications. In Section[6]we discuss future directions for the use of hierarchical data structures to visualise and 
analyse astronomy data. 


2. Hierarchical progressive surveys 

2.1. HEALPix 

HEALPix is a curvilinear partitioning of the sphere that supports a hierarchical tree structure for multi-resolution applications. The 
detailed geometry and properties of HEALPix are described in Gorski et al. (2005). Here we summarise the properties that are 
necessary for the definition of HiPS. 

The HEALPix partitioning of the sphere uses a base resolution that divides the sphere into 12 quadrilateral pixels, which are 
recursively divided into four self-similar smaller pixels in a 2x2 pattern to create finer and finer higher order meshes (see Figure[T|i. 
The resolution of the mesh at a given order k is defined by N S jd e , which is the number of divisions of the base resolution pixels and 
which doubles at each successive order, so that a HEALPix map of N S id e =2 k consists of \2N~ ide pixels. All of the pixels at a given 
order have equal area of £l pix — 7r/(3AF; w J, and the pixel centres are arranged in rings of equal latitude. HiPS uses the HEALPix 


1 p± _ '-1 _ sidzil _. 

nested numbering scheme (Gorski et al.|2005 ). 


2.2. HiPS tiles and pixels 

HiPS makes use of the hierarchical features of HEALPix to organise astronomical data into HEALPix maps of different orders. In 
this section we describe the HiPS tile and pixel structure for image data where both the tiles and their pixels are defined by HEALPix 
orders. This provides the general framework for HiPS that is later applied to other types of data. 

For image surveys we start by resampling the images onto a set of HEALPix maps of increasing order k. Each HEALPix map 
generated in this way is a representation of the original image data mapped onto a HEALPix mesh with a different pixel size, and 
each successive order has four times the number of pixels of the previous order. This multi-resolution representation of the original 
images provides the basis for visualising the data in a progressive way as the pixels that are required for a given view can be accessed 
from the pre-computed HEALPix maps, and the nested pixel numbering scheme provides a simple hierarchical indexing system that 
encodes pixel inheritance across the different orders. 

It would however be impractical to access and transfer individual HEALPix pixels, so the HiPS scheme groups pixels into tiles 
as the basic unit of a map that is transferred when any of its pixels are needed. The specification of the tiles makes use of the 
hierarchical properties of HEALPix to define the tiles as the HEALPix pixels of a lower order HEALPix grid, so that both the tiles 
and the pixels within the tiles are described by the HEALPix geometry. The general relationship between the tiles and pixels is that a 
tile with n t ne pixels along each side forms a HEALPix mesh of order of k t u e < k such that log 2 (n ft 7 e ) = log 2 (iV S M e /N S ide,tUe) = k-k t n e . 
The HiPS structure is therefore in general defined as a set of HEALPix maps with orders k and tiling The HiPS tiles and pixels 
are illustrated in detail in Figure [2] 

For astronomical purposes the generic HEALPix spherical coordinate system (6, cp as in Gorski et al.|2005 1 is mapped to astro¬ 
nomical coordinates so that the positions of the pixels have precisely defined positions on the sky (Figure[T|i. For HiPS we consider 
only two astronomical coordinate systems; the equatorial coordinate system defined by the international celestial reference system 
(ICRS), and also the Galactic coordinate system. Limiting the coordinate system to these two options is intended to provide a good 
level of consistency while also accounting for current use of HEALPix. While there has been significant use of HEALPix with the 
galactic coordinate system for cosmic microwave background (CMB) data, it is much more common for survey data in other wave¬ 
bands to be described by equatorial coordinates. ICRS is chosen as the default coordinate system for HiPS with Galactic coordinates 
as the alternative option. 

The HiPS representation of a collection of astronomical images is comprised of a set of HEALPix maps over a range of suc¬ 
cessive orders k mm to k max . The choice of the number of pixels in the tiles is based on the need for efficient transfer of the tiles 
to an interactive visualisation client. For HiPS that only use low orders with original data resolutions ^1°, this is not so important 
because the volume of data to be transferred is never very large, and these data are best served by tile sizes n,u e - 64 to 256. For 
higher resolution data, practical tests with bandwidth of ~1 Mbit s 1 shows that tiles of 512x512 pixels (i.e. n,n e = 512) are a good 
balance between the size of the files needed to store the arrays, and the number of pixels to be displayed and updated in the (typically 
1000x1000 screen pixels) window of a visualisation client tool. 


http://cds.unistra. fr 
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Fig. 2. HiPS tiles and pixels. A HiPS structure is defined by a set of tile orders k ti i e and by the number of pixels along each side of the tiles n nte . 
Here we show an orthographic projection of the sphere shaded to indicate different orders, and highlighting in black three individual HiPS tiles 
of orders k„/ c =3,4 and 5, with pixel tiling of n, W( ,=512. Note that HiPS representations of astronomical data usually start with L, ife =3 as the lowest 
order tile, and each successive order of smaller tiles is a higher resolution representation of the data. Some characteristics of the tiles are shown in 
the expanded view of the tiles on the right of the figure, including the approximate side length of the tile 6 ti[e (defined as the square root of its area 
6 ti i e = The number of pixels inside a tile is constant over all of the orders, and in this example where n,u e =512 each tile comprises 512x512 
pixels. The pixels inside a tile correspond to a higher order HEALPix grid, and in this example k plx =k,n e + logo n li i e =k,n e + 9, so that the k ti i e =3 tiles 
are filled with 51'.'55 pixels that correspond to a HEALPix grid of order 12. 


The properties of HiPS map tiling and the relationship of the tiles to HEALPix pixels are shown in Table |T| The first four 
columns of this table show the properties of HEALPix; k, N S ide, N P i X , and 0 p i X following Table 1 of Gorski et al. (2005) which has 


been extended here to show all orders up to k— 29. As described above, the HiPS tiles that enclose the pixels of a given HEALPix 
order are defined by a lower order HEALPix grid, and this is illustrated in the table with our commonly used case of n t u e = 512. The 
three next columns thus show; the HEALPix order of the tiles k t u e , 5 \ 2 , the number of such tiles that cover the sky and the 

size of these tiles 6,n e s\ 2 . The zeroth order grid for the n t u e = 512 HiPS (k,i / e ,5i2=0) is therefore comprised of 12 tiles of size ~58!6, 
where the 512x512 pixels of the tiles correspond to HEALPix pixels of order k-9. Note that the nine order difference between the 
HiPS tile grids and the HEALPix pixel grids means that the k t u e , 512 , N t u e , 5 i 2 , and fy, 7^512 columns are shifted versions of the k, N P i X , 
and 9 p i X columns. Use of different choices for n t u e would result in a different shift, for example for some low resolution data sets we 
sometimes employ »„/ e =64, for which the values of 64 , NuieM’ and QtUeM would correspond to a shift of log 2 n t u e = 6 rows in the 
table with respect to the k, N S ide, and N P j X columns. The properties of single order HiPS with n t u e = 64, 128, and 256 are shown in 
Table 2. 


2.3. HiPS file structure and indexing 

The tiles define the basic unit of storage for the set of HEALPix maps that make up a HiPS. Due to their quadrilateral nature, the 
HEALPix tiles are readily encoded as square arrays that can be stored in a number of different file formats. The simplification of 
having square arrays is a major advantage of HiPS as it leads to straight forward methods for drawing the tiles in any sky projection. 
It also means that the description of the array stored in the file is very simple in that all of the array positions are filled. 

The HEALPix indexing scheme provides a tree structured numbering system for HEALPix pixels with a simple relationship 
between the pixel numbers at different orders. Using the nested system leads to a simple organisation of the tiles into a filesystem 
directory structure. Each tile is encoded as a file, and all of the tiles for a given order are stored in a directory for that order. Starting 
at the third order (Aq, 7 ^ 2 = 3 ), the directories are named Norder3, Norder4... NordeiT,,/,,^ up to the maximum tile order. Within 
each Norder directory^ the tiles are named according to their tile number within that order i.e. NpixO, Npixl... Npix(A f , 7^512 - 1). 

The square arrays of pixel values in the tile files are arranged so that the pixels are in their natural spatial order. This means that 
the array when displayed as a square image, is a slightly distorted view of the curvilinear quadrilateral area of the sky. The natural 
spatial ordering of pixels in the tile means that no pixel indices or astrometry information needs to be stored in the file, and the 


simple ordering also facilitates the computation of the geometric projection of the pixels onto a display of the sky (see section 5.2 1 . 


10 The Npix files benefit from being grouped into further sub-directories of 10000 Npix files named DirlOOOO, Dir20000, etc. In this scheme, an 
individual pixel n at order k is in the tile m = n/( 2 10 ) at order k t n e , 512 = k-9, and will be found in the file Norder(£ - 9)/Dir(m, 10000)/Npixm. 
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Fig. 3. Multi-Order Coverage (MOC) map of the SCUBA 850 gm data. The lower left panel shows an approximately 5°x2?5 region of the HiPS 
image of the SCUBA 850um map. The MOC map of this region is shown in red, where we display the individual HEALPix pixels of different 
orders that define the MOC. The lowest order (largest HEALPix pixel) in this MOC is k= 6 (corresponding to B- 55;0), and the range of orders 
extends up to k=15 as shown in the magnified section of the figure where the smallest (0=6'.'44) pixels define the resolution of the MOC. 


Also, the use of the simple files and directories (compared to a database management system, or special format) means that HiPS 
are easily transportable and may be made available over the internet without the need for special servers or protocols. HiPS files 
may be simply put online with a standard web (HTTP) server. 


3. Multi-order coverage maps 

The HEALPix framework has also been used as the basis for the hierarchical multi-order description of sky regions called Multi- 
Order Coverage (MOC) maps. The MOC system exploits the same hierarchical characteristics of HEALPix as used by HiPS, in this 
case to describe regions of the sky in terms of coverage by HEALPix pixels of different orders. MOC has already been accepted as an 
IVOA recommendation (Fernique et al. 2014] > for the description of sky regions. The common underlying framework of HEALPix 
means that MOC and HiPS are both highly interoperable and complementary to each other. 

MOC maps can support the description of many types of sky coverage from a complete 4tt sr full sky survey, to the non¬ 
contiguous complex coverage of a set of pointed observations. The multi-order nature of MOC maps allows them to be defined at 
any HEALPix resolution to support different kinds of usage. Coverage maps that characterise the outlines of sky regions are very 
useful for many purposes such as establishing intersections, unions, and other logical operations based on the geometric coverage 
of regions of the sky. 

MOC maps are defined by first considering the set of HEALPix pixels at a given order that cover a region on the sky. This 
set of pixels is simply described as a list of HEALPix pixel indices in the nested numbering scheme. The self-similar hierarchical 
properties of HEALPix can then be used to compact the list into a concise and unique set of HEALPix pixel indices from multiple 
HEALPix pixel orders that provide the equivalent sky region. As such, MOC maps can be directly defined from any kind of data 
that is distributed on the sky. 

In practice a MOC map is constructed by selecting a maximum HEALPix order that will determine the size of the smallest 
HEALPix pixels to be used, and then listing all of the pixels that intersect the data. The list is then compressed by replacing any 
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four consecutive tiles at order k by their parent tile from order k — 1. This is done recursively down to k - 0. resulting in a list of 
HEALPix indices from different orders. The encoding of these lists of HEALPix indices into various formats including ASCII, and 
their serialisation into the FITS format, makes use of NUNIQ packing as described in Fernique et al. ( 2014| . 

A MOC description of a sky region is an approximation of the exact region boundaries. The resolution of a MOC is determined 
by the highest order used, and the outline is always limited to the compound shape of multi-order HEALPix pixels. Higher order 
MOC maps can be used to improve accuracy, but the limitations of the MOC approximation must be taken into account when using 
it for different purposes. While detailed schemes for expressing the exact coverage of sky regions are available, e.g. Rots ( |2011j >, 
the emphasis of MOC is on facilitating its ease of use and implementation. Another benefit of the simplicity of MOC is that the 
description of a sky region is unique, whereas schemes that use geometric constructions in various coordinates systems can be 
described in multiple (non-unique) ways especially in the case of complex regions. 

Figure [3] provides an example of a MOC that represents the coverage of a portion of the SCUBA (Submillimetre Common User 
Bolometer Array) 850/rm data (see P/SCUBA/850em in Table 3.). SCUBA is a bolometer camera on the James Clerk Maxwell 
telescope (JCMT) ( [Holland et al.|1999] >. A HiPS image of all of the public 850 pm observations has been constructed by the CADC. 
The combined data have an irregular pattern of coverage on the sky because of the distribution of the telescope pointings and the 
scanning characteristics of the array. The MOC representation of this data spans all of the regions where the data is defined, and 
provides a convenient footprint of the survey region on the sky. This MOC has been calculated by identifying all of the HEALPix 
pixels of order k = 9 (corresponding to a resolution of 6(871) which include valid data, and then compacting this list of pixel indices 
to create a unique list of HEALPix pixels over all orders k ^ 9. The resulting MOC displayed in Figure [5] shows how the HEALPix 
pixels of different orders fill out this complex and non-contiguous pattern of coverage on the sky. Compacting the pixels allows the 
larger contiguous areas to be efficiently filled with lower order pixels, while the finer details of the shape are represented by the 
smaller pixels at the chosen resolution. 


4. HiPS for images, source catalogues, and three-dimensional data cubes 

The development of the HiPS scheme has been driven by the need for scientifically robust hierarchical access to image survey data. 
The most direct use of HiPS is thus the mapping of astronomical images onto HiPS tiles. The HiPS scheme is however not limited 
to images, rather it can be used for many kinds of data. The concept that promotes the use of HiPS beyond its application to images 
is that the HiPS tiles may be used as general containers for any sort of information that is related to the sky coverage of the tiles. 
The content of the tile container can then vary according to different objectives. HiPS data structures have so far been successfully 
generated for astronomical source catalogues and multi-dimensional cube data, and the tile container has also been found to be very 
useful for storing links to the detailed metadata about the original progenitor data associated with the tiles. 

Here we describe the considerations to be taken into account when generating HiPS for different types of data; images, cata¬ 
logues, and three-dimensional data cubes. 


4.1. Images 

As described in Section [2] the HiPS representation of an image survey is generally constructed by resampling the images onto a 
HEALPix grid at the maximum desired order k max , and then generating the tile images for each of the tile orders. The process 
for resampling the pixel values from the original images depends on the intended purpose, and involves all of the same issues as 
generally encountered when mosaicking images. In general it is necessary to consider the desired angular resolution and resampling 
of the images, and also the choice of the methods to be used for combining data in regions where images overlap, and how to deal 
with variations in the background level. 

The choice of the highest order k max for a HiPS determines the minimum pixel size (see Table [TJ and is usually chosen to be 
close to the angular pixel size or resolution of the original data. To allow for various potential uses of HiPS no restrictions are made 
on the algorithms for resampling the original data onto the pixels, nor on the methods used for construction of the hierarchy of 
orders. The HiPS that are currently available for a wide range of surveys and image data sets demonstrate the benefits of some of 
the different approaches that have been tried so far. Many of these efforts have focused on the display quality of the final results, and 
others have placed greater emphasis on photometric accuracy and reducing re-sampling effects. These issues also guide the choice 
of file format to be used because of the different display efficiencies and dynamic ranges that may be supported by different file 
formats. 

The methods for generation of successive orders k max -1 to k mw of HiPS hierarchies may also be freely chosen depending on the 
intended purpose. The self-similar nature of the tiles means that conservation of flux in order k is achieved by simple averaging of 
each set of 2x2 sub pixels in order k— 1. The desired dynamic range of the pixel values influences the choice of the image file format. 
HiPS tiles are most appropriately stored as FITS format files when the full original dynamic range of the data must be preserved. 
For better performance however, JPEG or PNG encoding of the tiles is more appropriate. In this case with reduced dynamic range 
it is better to construct the successive HiPS orders using the median because it more reliably preserves the interesting structure in 
the images. For most cases we have chosen to generate both FITS and PNG format files. 

Construction of HiPS also requires consideration of the choice of the lowest order to include. Low order tiles (k t u e = 0,1, and 
2) correspond to large areas of the sky with tile side lengths 0 t n e > 10 deg. These orders are typically not useful for display of 
astronomical images because they are too coarse even for a full sky projection. Accordingly HiPS maps usually start at order k t u e =3 
and include all of the successive orders up to a maximum order chosen so that the HEALPix pixel size is close to the original image 
pixel size. 
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Our construction of a HiPS representation of the well known Digitized Sky Survevp] (DSS2) illustrates a number general 
aspects about creating HiPS for all-sky surveys. The digitized photgraphic plate images of the DSS2 (~ 1.5 million FITS format 
images of 768x768 pixels) have an original pixel scale of 170 pixel -1 and cover the entire sky. To map the DSS2 onto HiPS tiles 
of different orders we first identify the HEALPix pixel order that corresponds to the resolution of the survey, in this case £=18 with 
0pi X =0781 is a good choice as it slightly over-samples the 170 DSS2 pixels. A HEALPix map of the whole sky at this order contains 
^Vpu^S^SxlO 11 pixels. Choosing the HiPS tiles to be of our preferred size of 512x512 pixels means that the highest tile order will 
be £ri 7 e, 5 i 2 = 9 , corresponding to a tile size of 6, *,512=6(87. 

The most important step in the construction of the HiPS is the mapping of the DSS2 pixels onto the set of 0781 HEALPix pixels 
of order £=18. To re-sample the DSS2 pixel values onto the HEALPix pixels we use a bi-linear interpolation of the four closest 
DSS2 pixel values. When multiple overlapping input images contribute to a pixel we find that applying weights to input pixel fluxes 
as a function of the distance to the closest edge provides good results for the display quality of the tile images. (As the DSS2 pixel 
values are a measure of the photographic density of the original plate, which is non-linear with the intensity, we concentrate on 
display quality rather than photometric considerations in this example). Once the full set of £=18 pixels has been computed from 
the original DSS2 pixels it is straight forward to construct all of the lower orders by taking the median of each 2x2 pixels to create 
the successive orders. At each order we group the pixels into 512x512 tiles so that the tiles and pixels of the different orders have 
properties as shown in Table [I] The lowest tile order included is £ f *=3 (0,n e 512 = 7 : 33 ), which has pixels of HEALPix order £=12 
(0^=5175). 

The resulting HiPS of the DSS2 can be visualised on all scales from the whole sky down to the resolution of the survey by 
streaming the 512x512 pixel tiles of the required orders to a tool that can draw those tiles on a projection of the sky. Zooming and 
panning of the HiPS simply requires knowledge of the HiPS indexing structure so that the correct tiles for a given view can be 
accessed and mapped onto the display. 

A visually stunning example is the HiPS that has been generated for the GLIMPSE 36([^] survey of the plane of the Galaxy. 
GLIMPSE 360 is a part of the Galactic Legacy Infrared Midplane Survey Extraordinaire (GLIMPSE, |Benjamin et al.| ( |2003) >, 
Churchwell et al. ( |2009| ) observed by the Spitzer mission (Werner et a l.|2004 >. The GLIMPSE 360 survey comprises images at 3.6 
and 4.5 microns and the processed data are available from IRSAp^l A HiPS representation of the processed and mosaicked data has 
been generatecf 14 ] using the CDS Hips gen program. As the data were already stitched together as a high quality colour image, the 
construction of the HiPS concerned the mapping of the data onto the highest order grid of £=18, followed by computation of all 
of the orders down to £=12 with 512x512 tiles (corresponding to a range of tile orders k t u e =9 to 3). The result shown in Figure [4] 
illustrates how HiPS provides a whole-sky view of this data set, and how the tiles from different orders are used to visualise the data 
at different resolutions. 

For a number of the HiPS generated by the CDS there was a need to perform a background level adjustment of the original 
images. Simple subtraction of a constant background level has been applied to the 2MASS J-band survey using the original SKYVAL 
values. In other cases the background level has been estimated from a subset of pixels of each original image. This has been done 
for the SuperCOSMOS Her survey for which no background level was provided with the original data sets. Other examples required 
the use of masks to remove bad (or non-data) pixels from the original images. The GALEX AIS survey images have a circular data 
area within the rectangular data arrays (because of the circular GALEX field of view), requiring the use of masks to select the data 
values to generate the HiPS. 

For data sets that are already in HEALPix format, the generation of their HiPS representations is straight forward as it only 
requires derivation of the set of maps for the desired tile orders from the original HEALPix map. This is the case for maps of 
the cosmic microwave background (CMB) radiation that are available via the Legacy Archive for Microwave Background Data 
Analysis (LAMBDA^!- 

The WMAP Year 9 data, for example is provided by LAMBDA as a £=9 (A3;*=512) HEALPix map with 6 p j X = 6(87. This map 
can be directly expressed as a single order HiPS with £,*=3 and n t u e =64 by simply grouping the original HEALPix pixels into 
64x64 size tiles. Here there is no need for generation of lower orders because the data are already at the lowest useful HiPS tile 
order of three. 

The Planck LFI data have been provided by ESA/ESACj^jas £=10 (A3;*= 1024) HEALPix maps for the 20, 44, and 70 GHz 
bands, with the 70 GHz band also duplicated as a £=11 (N S id e = 2048) map. The HFI products at 100, 143, 217, 353, 545, and 857 
GHz are provided as £=11 (A3,*=2048) maps. As with WMAP, these Planck maps can be represented by single order HiPS with 
the £=10 LFI maps expressed as a £,*=3 grid with n,u e = 128. The £=11 LFI and HFI maps are naturally expressed as £,*=3 grids 
with n t u e = 256. See Table[2]for the properties of the HiPS tiles and pixels for HiPS with n t u e of 64, 128, and 256. 

Note that the original coordinate systems of all of the CMB data sets are Galactic coordinates. The HiPS of these data sets do 
not make any change in this respect as their HiPS representations simply use Galactic coordinate system HEALPix grids. 

The methods used to generate the HiPS described in this article are available in the CDS Hipsgen program. The preferred input 
files for Hipsgen are FITS images with valid WCS headers. Compressed (gzip, zip, hcompress) and multi-extension FITS files are 
also supported. The output formats can provide the complete dynamic range of the data or a more compact representation. The full 


11 The digitized sky surveys were produced at the Space Telescope Science Institute under U.S. Government grant NAG W-2166. The images of 
these surveys are based on photographic data obtained using the Oschin Schmidt Telescope on Palomar Mountain and the UK Schmidt Telescope. 
The plates were processed into the present compressed digital form with the permission of these institutions. See https://archive. stsci. 
edu/dss/faq.html 

http://www.spitzer.caltech.edu/glimpse360 

13 http://irsa.ipac.caltech.edu/data/SPITZER/GLIMPSE/ 

14 by Thomas Robitaille 

15 http://lambda.gsfc.nasa.gov 

16 European Space Agency / European Space Astronomy Centre 
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Fig. 4. HiPS of the Spitzer GLIMPSE 360 survey. The top panel shows a full sky Aitoff projection of the GLIMPSE 360 HiPS, rotated so that 
the galactic plane is horizontal and overlaid with a (green) grid of the lowest order k,n e =3 HiPS tiles with their tile index numbers. Note that only 
tiles that intersect the data are included in the HiPS and that regions of those tiles where there is no data are set to black. The lower panel shows a 
magnified view of the HiPS as shown by Aladin overlaid with the approximation of the HiPS tiles as they are projected onto the interactive display. 
The green grid shows how Aladin constructs the current view from k n7( ,=3 tiles, where in some parts of the view it is necessary to subdivide the 
tiles to gain more control points to minimise the distortion related to the bilinear tracing of the data onto the screen view. 


pixel values can be provided using FITS format. Using JPEG or PNG output formats provide a more compact result that degrades 
the dynamic range of the pixel values to 256 grey levels (or 255 grey levels plus one transparency channel). 

Colour images may be generated using two or three input data sets corresponding to RGB channels. It is possible to create 
coloured HiPS based on JPEG or PNG coloured tiles, either from an original coloured data set (for example the SDSS colour HiPS 
provided as a collection of JPEG images with an astrometric calibration), or from two or three already computed HiPS images 
where the individual tiles can be combined into colour RGB tiles. 

The volume of data in a HiPS representation of a full sky survey is typically an increase of ~ 1.5 times the volume of the original 
data. The total number of pixels over an entire HiPS structure will be limited to a factor of 4/3 (2^L 0 $r) greater than the /V ;)iA of 
the highest order map because of the four-fold decrease in the number of pixels within each lower order. For example, summing 
Npi X (from Table 1) over the k- 18 to 12 orders used for the DSS HiPS shows that the total number of HEALPix pixels in this HiPS 
is increased by a factor of 1.33325 over the k- 18 map. In practice the file storage of these pixels leads to a greater increase, for 
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example the 1.472 TB of DSS original files becomes 1.39 times larger as a 2.048 TB HiPS. The file formats have a much more 
significant effect on the volume, with JPEG and PNG formats being 20-40 times smaller than the FITS representation (dependent 
on image depth and structure). 

4.1.1. HiPS for collections of pointed observations 

The HiPS scheme is sufficiently general to handle data with any kind of sky coverage and there is no restriction on the contiguity 
of the input images. This facilitates the use of HiPS to combine collections of data, such as the pointed observations of a given 
instrument. The archives of the Hubble Space Telescope (HST) at the CADC have been used as a testing ground for the development 
of HiPS for pointed observations. Launched in 1990, the HST has acquired more than one million observations of which some 
887000 are direct images. Although HST is not a survey instrument, the quality and depth of most images are extremely valuable 
for a vast variety of astronomical projects, and we consider that the HiPS system for HST will open up new ways of using these 
images for research. 

To compute the HiPS for public HST images one has first to decide how those images should be grouped together. Consideration 
of the distribution of the observations in the many different filters (and filter combinations) shows that HiPS generated for a subset 
of ~15 filter groups would cover most of the HST observed fields. In defining these groups it makes sense to combine similar 
bandpass filters into a single group based on the number of images that use those filters and their characteristics. In addition, it is 
intended to eventually provide an achromatic HiPS for all the HST images. For each of the selected 15 biters groups we extracted 
the calibrated science data as input for the HiPS (using the CADC HST cache system). This was done for single images and for the 
already combined composite images (associations). 

Extensive trials were done to hnd the best parameters to generate this brst set of HST HiPS. The biggest challenges were to 
bnd the best choices to manage the sky brightness, and to handle the small belds of view. For sky brightness, we hnd that using 
an automatically estimated background value for each image produces a more uniform visual product. The HST instruments ACS, 
WFC3, NICMOS, and WFP2 have intrinsic angular resolutions between 50 mas and 071, so we tested different values of k max for 
the highest order maps of these data. Following these tests we chose to over-sample the data by using a maximum pixel order 
of k max = 23, which has 6 P i X = 25.1 mas (with corresponding tiles of order k,n e = 14, 0 t u e ,5n= 1271). Presently we use a single set of 
parameters to generate all of the HST HiPS, and this is done with the Hipsgen progranf* 7 ] 

During the trials we found that it was extremely useful to preserve the access links to the original images that contribute to a 
given HiPS tile, along with a small set of metadata that describes the images. Given the size of these data collections, it is important 
for the practical (data access) efficiency to extract this information while computing the HiPS. This metadata includes the exposure 
time, proposal number, observation date, and other information which is then stored with the tiles. This information can then easily 
be used alongside the tile images when visualising the HiPS, such as to provide a link (e.g. as an overlay) on each tile to the original 
data in the archive. 

Generating HiPS from large data collections such as HST requires significant planning and consideration of the growth and 
maintenance of the system. As the HST is still in operation, the CADC has decided to re-generate the HiPS for each of the filter 
groups every six months or so. Systems for incremental updates will likely be needed to streamline such operations. 

4.1.2. HiPS images and MOC 

MOC maps may be used to support HiPS data sets in a number of ways. A MOC map that has been constructed from a HiPS data 
set at the resolution of the smallest HiPS tile provides a convenient summary of the coverage of a HiPS data set. The construction 
of a MOC at any chosen HiPS tile resolution is straight forward because a HiPS map contains only the tiles that have data values, 
and the list of tiles may be used as the initial list of HEALPix pixels for the dehnition of its MOC. As described in Section 3, the 
initial list is then compacted to build the MOC. In this way, MOC maps of any chosen resolution may be constructed. 

We found that it was useful to generate a low-resolution MOC maps for each of the HiPS data sets hosted by the CDS, where the 
resolution of each MOC is chosen to be about four orders lower than the highest order (i.e. k ma t -4) of the HiPS. These MOC maps 
are stored with the metadata of the HiPS. The ease of generating such MOC maps has facilitated the implementation of interactive 
MOC generation tools in Aladin, which includes the capability to specify and change the MOC resolution. 

4.2. Catalogues 

The development of HiPS for catalogues is aimed at enabling the manipulation and visualisation of large (— 10 9 12 source) cata¬ 
logues. HiPS provides a means to organise a catalogue based on the spatial distribution of the source catalogue points on the sky. A 
HiPS catalogue made up of multiple tiles over a range of tile orders ( k t u e ) can provide different views of the catalogue that change 
as a function of the angular resolution. The freedom to choose which subset of the catalogue to include at any given order means 
that, in addition to the spatial organisation, the catalogue may also be organised over the different orders according to a hierarchy 
that can be based on any quantitative property of the catalogue, such as the source brightness, or held density, source redshift, 
or distance. The association of subsets of catalogue sources with HiPS tiles over a range of orders enables a progressive view of 
the catalogue so that the selection of sources that are displayed can change as a function of the zoom level. For example, a HiPS 
catalogue organised by source brightness can be presented so that the widest full sky view shows only the brightest sources, with 
fainter sources progressively appearing in the display as one zooms into smaller and smaller regions (Figure [5]). 

17 java -Xmx6000m -jar AladinBeta.jar -skyval=true -blank=0 -force -in=srcDir -out=trgDir -order=14 -img="reference_image" - 
fitskeys="INSTRUME OPT_ELE DETECTOR TIME_EXP PROP_ID PRODTYPE TARGNAME NMEMBERS OBS_DATE" 
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Fig. 5. HiPS catalogue tiles - content and visualisation. The content of a HiPS catalogue tile depends on how the source catalogue points are 
chosen to be distributed over the different HiPS orders. The upper panel shows three HiPS tiles of successive orders, with the lowest order tile on 
the left and sibling tiles to the right. The filled circles in the tiles represent catalogue sources that are ’stored' in each tile, and the upper set of tiles 
represents a sequential view of these orders. The lower panel shows the same three tiles but in this case we represent a cumulative view of the HiPS 
catalogue, where zooming into higher order tiles shows the sources in the tile plus the cumulative set of sources over all of the lower orders. 


In practice a HiPS catalogue is organised in the same way as a HiPS image survey, making use of the HiPS tiles to organise the 
spatial sky coverage of the catalogue so that the astronomical source positions in the catalogue are associated with the corresponding 
HiPS tiles. Whereas each tile of a HiPS image survey contains a 512x512 image, each tile of a HiPS catalogue contains a list of 
catalogue sources (with at least the RA and declination coordinates of the sources) in the form of catalogue rows that have been 
extracted from the full original catalogue. In the present implementation of HiPS catalogues the tiles are stored as ASCII tab 
separated value (TSV) files, and these tiles are organised into the file system directory structure of HiPS tile orders. The metadata 
about the columns of the catalogue are stored in a VOTable (Ochsenbein et al.|2011| document alongside the lowest order tile files. 


The construction of the catalogue tiles involves extracting rows from the original catalogue into the HiPS tiles over a defined 
range of orders. The main considerations when constructing the catalogue tiles are to decide on the approximate number of sources 
per tile, and the definition of the sorting key. The sorting key determines the distribution of the catalogue sources over the HiPS tiles 
as a function of the HiPS order, and the sorting key can be defined by any quantitative property. 


Given the sorting key, catalogue tiles can be constmcted by successively filling the tiles with sorted order sources up to an 
approximate limit for each tile, and then recursively moving to the next order to fill the four sub-tiles. Practical tests show that 
catalogue source densities of up to ~500 sources per tile allow for manageable zooming, but the tile source limit need not be 
constant, and in some cases it has been found to be useful to allow logarithmic or other non-linear increases in the number of 
sources per tile as a function of the order. We note that HEALPix aliasing patterns can occur when the original catalogue has 
heterogeneous distribution, and that this can be minimised by modifying the number of sources in the tiles based on the catalogue 
source density on wider scales. 


Filling the tiles in this manner means that each catalogue source is only included once in the HiPS catalogue. Visualising a HiPS 
catalogue by scanning through the orders provides a progressive view of the catalogue, and this can be controlled in a number of 
ways. For example, a sequential view of the tiles of different orders would show only the sources in each order. More commonly, a 
cumulative view shows all of the sources in the area of the tile up to a given order. This is illustrated symbolically in Figure [5] where 
the top panel represents a sequential view of three successive orders of a HiPS catalogue, showing the catalogue sources that are 
stored in each order. The lower panel shows a cumulative view of the same HiPS catalogue, where zooming into successive higher 
orders shows all of the sources in the area of the tile up to the current order. Note that construction of a cumulative view requires 
combining all of the sources from the current tile and all of its parent tiles. 


The 2MASS all-sky catalogue of point sources (Cutri et al .|2003a| contains ~470 million sources with measurements in the J, 
H, and Ks bands. A HiPS representation of this catalogue has been generated by mapping the positions of the sources onto HiPS 
tiles up to a maximum order of k t u e = 11, which corresponds to a tile size of H72. The sorting key for this HiPS catalogue is based 
on the infrared brightness of the sources where we have combined the fluxes in the J, H, and Ks bands. 


The largest HiPS catalogue that has so far been constructed is that of the Gaia universe model snapshot (Robin et al. |2012a) . 
This is based on GUMS-10 simulation of the expected content of the catalogue that will be generated from the ESA Gaia astrometric 
mission (Perryman et al. 2001). Gaia is expected to create a catalogue of ~10 9 stars with astrometric accuracies of 5-15 //as (for 
G band magnitudes <12) and 5-600 //as for the full catalogue. The ~2xl0 9 simulated Galactic objects in GUMS-10 catalogue of 
"milky way stars" are based on the Besan£on model (Robin et a l,|2003| , which provides the distribution of the stars, their intrinsic 
parameters, and their motions. The HiPS of this catalogue uses the simulated sky coordinates and the sources are distributed over 
HiPS orders up to k fI / e =ll. The sorting key is chosen to be the simulated barycentric distance (in parsecs) that is provided in the 
catalogue, a quantity that the mission is expected to derive from accurate parallax measurements. 
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t, X, v, E,... 

Fig. 6. FiiPS cubes. FliPS cubes are simply a set of nested FliPS images, where the third dimension of the cube may represent time (t), wavelength 
(d), frequency (v), Energy (£), or any other quantitative parameter. HiPS stores the third dimension planes in the same way as FliPS tiles, using 
simple files, where the successive planes are labelled with a sequential suffix number. 


Another example is the HiPS catalogue that has been constructed from the CDS SIMBAEp^] astronomical database (Wenger 


et al.] [2'000| . SIMBAD provides basic data, cross-identifications, bibliography, and measurements for astronomical objects. SIM 
BAD currently includes over 7.6 million well documented astronomical objects and the database can be queried by object name, 
coordinates, and various criteria. While SIMBAD is not a single catalogue, much of the information can be expressed in the form 
of a catalogue of sources so that SIMBAD may be represented by a HiPS catalogue. Experimenting with different sorting keys for 
distributing the SIMBAD sources over a range of HiPS orders showed that sorting by properties such as source brightness, spectral 
properties, angular size, or literature statistics all provide new and useful ways of viewing SIMBAD. The publicly available SIM¬ 
BAD HiPS catalogue is sorted by the number of literature references associated with each astronomical object. The resulting HiPS 
employs tile orders up to k,u e = 14, which has a tile size of 1279. When viewed in Aladin the SIMBAD HiPS shows the most highly 
referenced famous sources in the zoomed-out wide scale view, with more sources appearing as the display is zoomed into higher 
orders. 

The methods used for generating the HiPS catalogues described above are all included in the CDS Hipsgen-cat program. 
Hipsgen-cat accepts inputs in the form of list of sources in CSV, FITS, or VOTable format, and has been successfully tested 
with catalogues of up to 2 billion sources. The typical generation time is a few hours for catalogues of the order a hundred million 
rows, for a maximum HiPS tile order of ~ 11. 


4.3. Cubes 

The development of HiPS cubes extends the concept of HiPS images into a third dimension to enable the hierarchical organisation 
of three-dimensional data cubes which have two spherical coordinate axes plus a spectral (wavelength, frequency, energy), velocity, 
or temporal axis as is symbolically illustrated in Figure[6] The basic idea here is that each channel (or plane) of the third dimension 
of the data cube is stored in additional HiPS tiles. 

The hierarchical description of the spatial part of HiPS cubes is identical to that of HiPS images. Experimentation with different 
axis compression schemes for the third dimension using a variety of different cube data showed that it is preferable to retain all of 
the third dimension planes, even when the spatial dimensions are heavily smoothed as occurs for low order HiPS tiles. This choice 
is mostly guided by tests of spectral cubes, where we found that visualisation of the sharp spectral features in the cubes is important 
even when the spatial dimensions were compressed, and that the blurring of spectral features in compressed data negated the benefits 
of the zoomed-out views of the cube. 

The number of third dimension planes in astronomical data cubes can vary greatly; cubes constructed from multi-band images 
may have only a few planes, whereas radio spectral data cubes and integral field spectroscopy cubes may have hundreds or thousands 
of planes. To avoid streaming of large volumes of data from HiPS cubes, as might occur for the low orders of HiPS cubes with many 
planes, we find that sampling of the third dimension to access only every nth (2nd, 5th etc.) plane is a satisfactory solution for low 
HiPS order projections. As such the control over the effective reduction in the number of planes to be downloaded is given to the 
client programs. 

18 http://simbad.unistra.fr 
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This approach also effectively prioritises fast spatial manipulation of the data ahead of spectral operations that might be per¬ 
formed by orthogonal cuts or extractions along the third dimension. By unpacking the cubes into separate frames we enable fast 
access to the individual planes of the cube. We have tested an alternative storage of the cubes in more integrated cube formats such 
as MPEG, which would allow for faster access (and compression) along the third dimension, but we found that the extraction of 
individual frames from multiple MPEG formatted cubes of adjoining tiles was too slow for the spatial manipulation of the cube. 

The construction of HiPS cubes is similar to that of HiPS images, with the extra dimension taken into account via the use of 
multiple frames for each tile image. These frames are simply multiple tile images in separate files and enumerated using a sequential 
suffix number, 0 to N(cube planes)-1. As is the case for HiPS images, the cube frame tiles may be stored as FITS, JPEG, or PNG 
files. For data that is already in cube format we unpack the planes to produce HiPS tiles for each plane. For multi-band surveys 
HiPS generation is straight forward, but care must of course be taken with the astrometric correspondence of all the band images. 
As such, HiPS of multi-band surveys may be simply constructed a posteriori from the HiPS of the individual bands. 

A HiPS cube has been generated for the Canadian Galactic plane survey (CGPS, Taylor et al. ( 2003| l). The CGPS is a wide 
area (1600 square degrees) atomic hydrogen survey covering a large fraction of the Galactic plane with an angular resolution as 
small as T. The original observations comprise 373 fields with the Dominion Radio Astrophysical Observatory (DRAO) synthesis 
telescope, and the processed calibrated daUrjconsist of sixty-eight 5°x 5°(1024xl024 pixel) data cubes with 272 spectral channels 
with velocity resolution of 1.3 km s -1 . The HiPS cube of this data was generated by mapping all of the image planes of the data 
cubes onto 272 HiPS structures using tiles of size 512x512, with a range of orders up to k max =14 corresponding to a maximum 
resolution of 1279. The highest order tiles are k t n e =5 which have an approximate size of G t n e , 5i2=l?83. The 272 HiPS structures are 
stored in an integrated way by combining all the files into a single HiPS structure and enumerating the tiles (with suffixes 0-271). 
The result is a homogeneous wide area HiPS cube that can be displayed all at once in a tool such as Aladin. The HiPS contains 
precomputed image tiles of all of the CGPS cube planes over the range of orders, so in addition to the wide scale view, we may 
zoom into the cube to access the higher resolution tiles, and at any order we may also interactively scan through all 272 spectral 
channels. In this way HiPS enables the CGPS data set to be accessed and visualised more easily than by loading the sixty-eight 
fixed-resolution individual data cubes. 

HiPS cubes have been constructed for the Calar Alto legacy integral field area (CALIFA) survey ( [Sanchez et al.|2012j >. The CAL- 
IFA survey is an on-going integral field spectroscopy survey of 600 local (0.005< z <0.03) galaxies obtained with the PMAS/PPak 
integral field spectrograph mounted on the 3.5 m telescope at the Calar Alto observatory. A HiPS cube was built from the 100 
low spectral resolution (V500, R-850) cubes of the first data release (DR1, Husemann et ak] |2013|l , and another HiPS cube was 
built from the extended (and superseding) sample of 200 galaxies in the second data release (DR2, Garcfa-Benito et al.|2014[ ). The 
CALIFA observations have 74"x64" hexagonal fields-of-view with 331 science fibres that project 277 diameters on the sky with 
gaps between the fibres that result in a filling factor of 0.6. The processed data cubes have l"xl"angular sampling, and a linearly 


sampled wavelength axis with 1877 channels that span the wavelength range 3749-7501A. The HiPS cubes are constructed from 
the processed data, taking advantage of the consistent format and cube dimensions to generate HiPS images for each of the 1877 
channels, and combining these into a single HiPS structure. Tiles of 512x512 pixels are used, and the maximum order of the tiles 
was chosen to be k,u e =9 corresponding to a resolution of 0781 to slightly over-sample the processed data. In the few cases where 
multiple CALIFA cubes overlap on the sky the average flux values are used. 

The CALIFA cubes cover only a small fraction of the sky, but their combination into single whole-sky HiPS structures (one for 
each of DR1 and DR2) brings about a number of benefits. HiPS provides a convenient way to access all of the CALIFA data with 
one action to reveal the distribution of CALIFA observations across the sky. The supporting metadata, including the sky coordinates 
of the original data cubes, provides a way to easily select and zoom into the sky regions covered by the data. Zooming into the 
higher orders of the HiPS cubes shows the details of the data including the spectral channels. Furthermore, multiple views of the 
same HiPS structure allows direct comparison of the different CALIFA galaxies, and allows one to interactively scan through the 
spectral channels of multiple cubes at the same time. In addition to the positions, the metadata also includes the links to the original 
data so that the full individual cubes may be loaded alongside their HiPS representations. Also the MOC that is included in the 
metadata can be used to make general queries for other data that intersects the same region of the sky. 

HiPS cubes have been constructed for a number of multi-band image surveys, where the planes of the cubes are the wavelength 
ordered bands of the survey. A good example is the HiPS cube of the IRIS (improved reprocessing of the IRAS survey) data 
( jMiville-Deschenes & Lagache|2005] l. To construct this HiPS cube the HiPS images of the individual 12, 25, 60, and 100 gm bands 
were combined into a single HiPS cube structure. 


5. HiPS data sets and tools 

5.1. Currently available HiPS data sets 

There are presently -200 publicly available HiPS data sets. Lists of the current HiPS images, HiPS catalogues, and HiPS cubes are 
provided in Tables 3, 4, and 5 respectively. 

Table 3 represents -180 HiPS image data sets. These are organised in the table by waveband: gamma-ray; X-ray; ultraviolet; 
optical; infrared; radio bands; and categories for emission line surveys of Hu, H i, and CO. Each HiPS has a unique CDS identifier. 
The table shows the base names of these identifiers, and each row represents either a single HiPS or a group of related HiPS. For 
example, the first row describes the four HiPS images of the Fermi data, one for each of the three energy ranges (3-300 MeV, 
300-1000 MeV, and 1-3 GeV) and another one for the colour composite HiPS of this data. 
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Table 3 lists the properties of each HiPS: the maximum tile order k,u e \ the size of the smallest pixel <t plx ; and the coverage of the 
HiPS as a percentage of the whole sky. Note that each HiPS image structure contains all of the tile orders from the maximum order 
down to k,n e =3 , so for example, each of the XMM HiPS with maximum values of k,u e -l are comprised of 5 successive tile orders 
&ri/e=3,4,5,6, and 7. The deepest HiPS in terms of the number of orders is that of the HST archive images, which has a maximum 
tile order of 14 and corresponding pixel size of 6 p i X = 25.1 mas. 

HiPS for which the maximum tile order is k,u e =3 are single order HiPS of relatively low angular resolution (~l'-7') such as 
the gamma-ray surveys, early X-ray surveys, CMB maps, and some radio astronomy surveys. There are 71 such single order HiPS 
(within 19 groups) in Table 3. 

There are 77 full-sky HiPS in Table 3. The highest order HiPS that cover the whole sky are the 2MASS and DSS HiPS with 
maximum tile order of k t n e =9 (N s ,d e -2 18 ) with 0 ; , il =O781. There are 55 HiPS that cover less than 5% of the sky, and 7 HiPS with 
less than 1% sky coverage. 

Table 4 lists the HiPS catalogues that are currently available. All of these have been generated from catalogues and data from 
the CDS Vizier and SIMBAD services. For each of these data we tabulate the CDS identifier, a brief description, and references to 
the original published catalogues. 

Table 5 lists the data sets that have been converted into HiPS cubes. These represent a wide range of cube characteristics in 
terms of the number of planes, field of view size, and distribution over the sky. 

Most of the HiPS in Tables 3, 4, and 5 have been created by the CDS from original or processed data that have been obtained 
from publicly available archives, or via collaborations with individual authors, projects, or partner data centres. There are 35 HiPS 
that are hosted outside the CDS, and there are 20 HiPS that have been computed externally to the CDS but that are hosted by the 
CDS. The table notes specify the sources of the original data and identify the parties who created the HiPS representations. We 
also also indicate the data centres where each HiPS is hosted. References to the published articles on the surveys, instruments, data 
releases, data processing, and data characteristics are also provided where possible. 

The lists of the currently known HiPS images, catalogues, and cubes are maintained at the CDS in the HiPS director)^] The 
HiPS directory provides information on each of the HiPS similar to what is presented in Tables 3,4, and 5, along with other practical 
information including the full CDS identifier and the root URL of each HiPS. The file formats that are available for each of the 
HiPS are also indicated, along with the date of last modification. The HiPS directory thus provides an overview of the current status 
of all HiPS data sets. 


5.2. Tools for accessing and visualising HiPS 


There are currently three astronomical visualisation tools able to visualise HiPS images. These are the Aladin desktop Java appli¬ 
cation and the Aladin Lite Javascript web widget both developed by the CDS, and Mizar which is a WebGL client developed by 
CNES as part of their SiTools 2 publication framework ( Malapert & Marseille|2011) . 

To visualise HiPS images a tool needs to perform a number of steps. Starting from the minimal information of the HiPS root 
URL, a tool needs to read the properties file of the HiPS to retrieve metadata about the maximum order, the size of the tiles, and 
the coordinate frame (ICRS or Galactic). To display the HiPS for a region of the sky the tool must compute the tile order that is 
best suited to the desired field of view. Then, using the HEALPix library ( |G6rski & Hivon||2011] >, and knowledge of HiPS tiling, 
the tool needs to identify which tiles it needs for the view, and to download the files for those tiles. To draw the tiles on the display, 
the display coordinates of the four corners of each tile must be determined, and then the tile image can be distorted through affine 
transformations to map it onto the display area defined by the four mapped points. 

Aladin is a multi-featured Java application for the discovery, access, and visualisation of astronomical data with detailed func¬ 
tionality for images and catalogues. Aladin is VO compatible, and can serve as a portal to access data made available by VO 
protocols. The early development of Aladin is described in Bonnarel et al. ( |2000) and current information is maintained at the 
Aladin web page^j The HiPS visualisation capabilities of Aladin have been developed alongside the development of the HiPS 
system at the CDS, and Aladin is the primary tool for testing and verifying HiPS. Aladin provides an interface to access all of 
the HiPS described in the HiPS directory, and it can also be used to access any available HiPS given its local file address or URL 
address. Aladin provides a familiar all-sky zooming and panning interface that makes full use of the multi-resolution features of 
HiPS. Multiple HiPS may be viewed simultaneously in separate windows, or layered together with transparency overlays, and views 
may be locked together for synchronised zooming and panning. For HiPS cubes, Aladin provides a facility for stepping through the 
third dimension of the cube plane by plane or as an animation. HiPS catalogues may be overlaid on images and zoomed to view 
the different HiPS catalogue orders. The general Aladin catalogue overlay feature also provides a very useful mechanism for going 
from a view of the HiPS image, to the original data sets that are listed in the progenitor metadata catalogues that can be included 
in HiPS. Aladin also provides a set of functions for using MOC maps (e.g. Figure 3 shows screen shots of Aladin being used to 
overlay a MOC map on the HiPS of the SCUBA data) including a multi-resolution display and the ability to combine MOC maps 
with logical operators to form unions and intersections. 

In addition to its capabilities for using HiPS and MOC, Aladin also includes all of the tools for creating HiPS and MOC maps 
from images, catalogues, and cube data sets. All of the HiPS described in this paper have been generated with the Hipsgen tool that 
is provided in the standard Aladin desktop application. Hipsgen can be used from the command-line or directly from the Aladin 
tools menu. 

Aladin Lite is a simplified HiPS visualiser that is developed in the JavaScript programming language. It runs in a web browser 
(without the need for a plug-in) and is able to display HiPS images (in JPEG or PNG formats), and HiPS catalogues with zooming 
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and panning capabilities similar to the Aladin desktop application. Other tabular data and information can also be overlaid. Aladin 
Lite can be used to visualise local HiPS data and HiPS that are hosted on any web server including those in the HiPS Directory. 
A strength of Aladin Lite is that it can be embedded easily on a third-party web page (Boch & Ferniqu e|2014j > and it comes with 
an application program interface (API) that provides the necessary controls to customise Aladin Lite to different needs. A good 
example of this is the Spitzer web pag^Jof the GLIMPSE360 project where Aladin Lite is used (alongside WWT) to visualise 
a colour HiPS of the GLIMPSE360 survey data. Other examples include the Aladin Lite implementations at CADlrjto preview 
HEALPix maps of various surveys, and the ADSp^all sky surve\p^| 


6. Discussion 

Concepts for the hierarchical organisation of astronomical data are an important part of how astronomy approaches the challenges 
of big data. The main ideas are that i.) hierarchical systems describe the data in a multi-resolution manner that facilitates efficient 
access to the data, ii.) hierarchical systems scale easily so that they may handle large and small data sets in the same framework, 
and iii.) the use of a generic hierarchical system for the wide variety of astronomy data enables new levels of interoperability. 

HiPS has been built with these underlying concepts in mind, and the result is a practical multi-resolution HEALPix data structure 
for astronomical images, catalogues, and three-dimensional data cubes. Here we discuss various visualisation and scientific aspects 
of HiPS. An important distinction of HiPS is that it goes beyond visualisation to take the scientific considerations of astronomical 
data into account providing a new statistical view of the data. Also, HiPS is designed to be easy to implement and we outline how 
it can be used to enable data sharing and interoperability. 

6.1. Scientific visualisation 

The description of HiPS in this paper has been largely presented from the point of view of visualising astronomical data. This is a 
convenient way to describe HiPS and indeed a number of the characteristics of HiPS, such as the preferred tile size, are motivated 
by practical considerations related to accessing and streaming of multi-resolution astronomy data to a visualisation tool. In terms of 
supporting visualisation of astronomical images, HiPS has similar performance to the other hierarchical formats that are used for this 
purpose (TOAST/HTM). The HiPS system based on HEALPix, provides a means to conserve the scientific properties of astronomy 
data, and goes beyond the description of images to extend the hierarchical description of data to catalogues and three-dimensional 
data cubes so that these data types may also be visualised as multi-resolution structures. 


6.2. HiPS facilitates scientific use of astronomical data 


To make HiPS a hierarchical framework that supports scientific use of data we considered that it was important to i.) use a sky 
tessellation that was directly translatable into astronomical coordinates, ii.) support accurate photometric properties of images that 
would allow scientific computation directly on the hierarchical grid pixels, and iii.) provide a way to link directly to the original 
progenitor pixel or catalogue data. 

HEALPix has proved to be a good choice because it automatically supports the first two considerations, and the catalogue 
and metadata support mechanisms of HiPS provide a simple way to include links to the original data. Also, HiPS benefits from 
the scientific computing libraries and tools that have already been developed for HEALPix. These libraries and tools are largely 
from the CMB and extended-source astronomy community, but there is also increasing use of HEALPix for other projects, for 
example Gaia is using HEALPix (private communication William O’Mullane; Mignard|2013j l. There is much potential for the wide 
spread use of HiPS and integration of HiPS functions in astronomy software. This would build on our earlier work to define MOC 
(Fernique et ak] [2014| that has resulted in MOC functions being included in the official HEALPix library ( jGorski & Hivon]201 1[ 
private communication Martin Reinecke). 

HEALPix, and hence HiPS, provides a direct translation from the pixels to the astronomical coordinates. HiPS supports both 
Galactic and equatorial (ICRS) coordinate systems. This duality is in response to the fact that many data sets have a natural coor¬ 
dinate system of Galactic or equatorial, leading us to support both systems. All of the current HiPS visualisation tools support both 
systems so that Galactic and equatorial HiPS can be used together. Computational interoperability between HiPS however requires 
the conversion of the data to one system or the other, and there may be cases where it will be necessary to host both Galactic and 
equatorial versions of HiPS. 

HiPS manages both visualisation and scientific computational considerations by allowing the HiPS tiles to be encoded in both 
visualisation and science data formats. HiPS tiles can be duplicated in multiple formats, and these are simply stored alongside each 
other in the HiPS file structure. FITS files are typically used for science grade data with the full original dynamic range of the data 
values encoded in FITS in the standard manner. The efficient transfer and display features of the JPEG and PNG formats facilitate 
their use in visualisation systems. Having both means that one can visualise the data efficiently, and then seamlessly switch to 
the strict science data format to access the actual measurement values and perform computations. Other formats or compression 
schemes (e.g. McEwen et a l.|201 1| ) could be explored as alternative tile encodings. 

Another way in which HiPS facilitates the interaction between data visualisation and analysis is by providing a mechanism for 
direct links to the progenitor original data. The information about the original data is best included into the HiPS at the time of 
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generation (especially for large data sets). This facilitates making catalogue style overlays of clickable links to the original archived 
data. This is also supported by tools that can display the original FITS data alongside of the HiPS representation. 

In terms of the photometric accuracy of HiPS, the successive orders generated by simple averaging (as usually done for FITS 
format HiPS) are strictly flux conserving. The initial construction of the highest order map is of course a re-sampling of the original 
data (except for data that is originally in HEALPix format), and the photometric characteristics of the data are determined by that re¬ 
sampling. While the tools provided in Aladin for generating HiPS ( Hipsgen ) provide basic resampling and mosaicking capabilities, 


there can be advantages to pre-mosaicking data with more control over the photometric properties using tools such as SWarp (Bertin 
|et al.|2002 1 or Montage (Jacob et a l .|2Q 1 Q[ >. We also note that Paradis et al. (2012) have already made progress with drizzle techniques 
for HEALPix data. 


6.3. Interoperability of images, cubes, catalogues, and coverage maps 

HiPS and MOC allow a new level of interoperability between image surveys, data cubes, catalogues, and sky coverage maps. The 
common basis of HEALPix means that both large and small scale data sets may be easily visualised and compared. The expression 
of sky regions as MOC maps facilitates the comparison of sky regions covered by different data sets (of any type) such as to establish 
regions of intersection between multiple surveys. With MOC, such comparisons reduce to the simple comparison of lists of integer 
HiPS tile numbers, and this method scales very well when there are many data sets to be compared. MOC is also very useful 
for establishing whether a given sky coordinate falls within the coverage of a survey, or within the instrument fields of view of 
observations in an archive. These capabilities also extend to the computation of the intersections between survey areas and sources 
in large catalogues. 

Furthermore, services that provide access to data with a given sky coverage may employ a MOC description to vet in-coming 
(coordinate based) queries. A service may make an initial comparison of the coordinates of an in-coming query with the MOC 
map coverage of the data. This allows the service to provide a fast initial response that would indicate whether a detailed database 
query of the service is likely to yield a result, hence avoiding complex queries on services that would return null results because of 
un-matched sky coverage. 

6.4. Implementation of HiPS 

The use of a simple file and file system structure to organise the tiles of a HiPS makes it very easy to publish HiPS to the web. The 
simplicity of such systems is a very practical consideration as it significantly lowers the barrier to implementing HiPS compared to 
systems where special databases or web servers are required. The use of simple files also means that the system is transportable and 
scalable, and can be easily duplicated on multiple servers (for example as mirror copies). Such transportability facilitates the local 
generation and testing of HiPS before deployment on public systems. HiPS is scalable in that the individual component files are 
small (i.e. managing large HiPS data sets involves managing large sets of small files) so that increases in volume are not affected by 
file size limitations. 

With the growing number of HiPS currently being generated by different data centres, projects, and individual astronomers, it 
is important to be able to organise and find these HiPS. Presently CDS maintains the HiPS directory, but there is clearly a role for 
virtual observatory registries to enhance the description and discovery of HiPS data sets. 


7. Conclusions 

We have introduced HiPS as a practical multi-resolution HEALPix data structure for astronomical images, catalogues, and three- 
dimensional data cubes. There are already many (— 180) HiPS image data sets available for a diverse range of sky survey data over 
all wavebands, and with sky coverages that range from a few square degrees up to the full sky, demonstrating the applicability this 
hierarchical method to many areas of astronomy. 

The development of HiPS has placed a strong emphasis on conserving the scientific properties of the data alongside both 
visualisation considerations and emphasis on the ease of implementation. The scientific robustness of HEALPix enables HiPS to be 
a framework for multi-resolution statistical analysis of astronomical data sets in that each level of a HiPS is a statistical summary 
of the information in the lower levels. Also, there is much potential for the development of astronomical analysis procedures that 
operate directly on the HiPS pixel and tile grids, benefiting from the computational libraries developed for all sky HEALPix maps, 
that may lead to new capabilities for analysing partial sky maps. In addition to the scientific benefits related to HEALPix, HiPS also 
provides a simple mechanism for directly linking to the original data. 

The extension of HiPS to three-dimensional data cubes provides methods to visualise large data cubes and to manage many 
cube data sets simultaneously. The HiPS cube examples presented here illustrate the use of HiPS for cubes with a spectral axis as 
the third dimension, but this dimension may also be used to represent time or other kinds of axes. 

HiPS catalogues use the HiPS tile structure to organise astronomical source catalogues. We have shown that this method is a 
practical way of handling very large catalogues with our examples demonstrating good performance for catalogues with up to 2xl0 9 
sources. The development of HiPS catalogues is at a relatively early phase with the eleven examples presented here representing our 
experimentation with catalogues of different size, density, and for which we have tested different sorting keys. These techniques will 
be developed further as they are expected to be very important for future large catalogues that will be produced by GAIA, LOFAR, 
LSST, Euclid, and SKA and its pathfinders. 
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The combination of HiPS and MOC provides a very powerful tool for astronomy because they enable new levels of interoper¬ 
ability between astronomical images, data cubes, catalogues, and the description of regions of the sky. The simplicity and ease of 
use of the HiPS system is designed to facilitate wide spread implementation of HiPS across all areas of astronomy. 
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Table 1. HiPS tile and pixel properties for n„7 e =512 


k 

N side = 2* 

N pix 

@pix ktile, 512 

Ntile, 512 

0tile,5n 

0 

1 

12 

58?6 




1 

2 

48 

29?3 




2 

4 

192 

U°.l 




3 

8 

768 

7?33 




4 

16 

3072 

3?66 




5 

32 

12,288 

1?83 




6 

64 

49,152 

55:0 




7 

128 

196,608 

27i5 




8 

256 

786,432 

13:7 




9 

512 

3,145,728 

6:87 

0 

12 

58?6 

10 

1024 

12,582,912 

3:44 

1 

48 

29?3 

11 

2048 

50,331,648 

i:72 

2 

192 

14?7 

12 

4096 

201,326,592 

5175 

3 

768 

7?33 

13 

8192 

805,306,368 

2578 

4 

3072 

3?66 

14 

2 14 

3.22 xlO 9 

1279 

5 

12288 

1783 

15 

2 15 

1.29 xlO 10 

6744 

6 

49152 

55:0 

16 

2 16 

5.15 xlO 10 

3722 

7 

196608 

27:5 

17 

2 17 

2.06 xlO 11 

1761 

8 

786432 

13:7 

18 

2 18 

8.25 xlO 11 

0781 

9 

3,145,728 

6:87 

19 

2 19 

3.30 xlO 12 

0740 

10 

12,582,912 

3:44 

20 

2~o 

1.32 xlO 13 

0720 

11 

50,331,648 

i:72 

21 

2 21 

5.28 xlO 13 

0710 

12 

201,326,592 

5175 

22 

2 22 

2.11 xlO 14 

50.3 mas 

13 

805,306,368 

2578 

23 

2 23 

8.44 xlO 14 

25.1 mas 

14 

3.22 xlO 9 

1279 

24 

2 24 

3.38 xlO 15 

12.6 mas 

15 

1.29 xlO 10 

6744 

25 

2 25 

1.35 xlO 16 

6.29 mas 

16 

5.15 xlO 10 

3722 

26 

2 26 

5.40 xlO 16 

3.15 mas 

17 

2.06 xlO 11 

1761 

27 

2 27 

2.16 xlO 17 

1.57 mas 

18 

8.25 xlO 11 

0781 

28 

2 28 

8.65 xlO 17 

0.786 mas 

19 

3.30 xlO 12 

0740 

29 

2 29 

3.46 xlO 18 

0.393 mas 

20 

1.32 xlO 13 

0720 


Table 2. HiPS tile and pixel properties of single order HiPS with n,u e =64, 128, and 256 


Utile-64 

k 

Nside 

Npix 

Opix 

ktile, 64 

Ntile, 64 

@tile,64 

9 

512 

3,145,728 

6:87 

3 

768 

7733 

Utile- 128 

k 

Nside 

N P ix 

Opix 

ktile, 128 

OO 

<N 

,$r 

Otil e ,m 

10 

1024 

12,582,912 

3:44 

3 

768 

7733 

futile — 256 

k 

Nside 

Npix 

Opix 

ktile,256 

Ntile,256 

Otile,256 

li 

2048 

50,331,648 

Y.12 

3 

768 

7733 
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Table 3. HiPS images 


ID" 

Data set name and description 

ktile 

Qpix 

Sky 

(%) 


- Gamma Ray - 




P/Fermi £ bcd . 

Fermi* 11 : 4 HiPS in 3-300 MeV, 300-1000 MeV, 

3 

6:87 

100 % 


and 1-3 GeV and a colour composite (n t u e = 64) 




P/EGRET £bcd . 

EGRET (2) : 12 HiPS for 100 MeV inf and sup maps, and 

3 

5175 

100 % 


diffuse maps 30-50, 50-70, 70-100, 100-150, 150-300, 





300-500, 500-1000 MeV, 1-2 GeV, 2-4 GeV, 4-10 GeV 





— X-ray - 




p/xmm £ ^. 

XMM-Newton : 5 HiPS for EPIC stacked images, PN 

7 

3722 

7% 


in 0.5-1, 1-2, and 2-4.5 keV bands and PN colour 




P/INTEGRAL £ s crf . 

INTEGRAL* 31 SPI GC (20-40 keV) 

3 

5175 

17% 

p/intgal £ cdh . 

INTEGRAL IBIS (4) : 4 HiPS for 17-35, 17-60, 

3 

5175 

31% 


35-80 keV, and colour composition 




P/RASS 0 c d ' . 

ROSAT X-Ray All Sky Survey* 51 

4 

2578 

100 % 


- Ultraviolet - 




P/GALEXGR6 f; ' rf/ . 

GALEX* 61 : Colour composition of GR6 AIS 

8 

1761 

69% 


- Optical - 




P/DSS2 £e ' i ' / . 

DSS2: 2 HiPS for red band and colour composite 

9 

0781 

100 % 

P/SDSS9 £ec,ir . 

SDSS DR9* 71 : 3 HiPS for r, g, and a colour composite 

10 

074 

36% 

P/Mellinger/color £c "' m .. 

Mellinger optical survey 

4 

25777 

100 % 

P/CFHTLS/D £ cd " . 

CFHTLS* 81 Deep : 6 HiPS for u,g,r,i,z, and colour 

12 

071 

0 .01% 

p/cfhtls/w £ cdn . 

CFHTLS* 81 Wide : 6 HiPS for u,g,r,i,z, and colour 

11 

072 

0.4% 

P/HST £ ^« . 

HST: 17 HiPS forNICMOS, WFPC2, and ACS in 

14 

25.1 mas 

0.044% 


FI 10W, F160W, F225W, F300W, F450W, F606W, 





F625W, F702W, F814W, F850LP, and GOODS* 91 





associations in b,v,i,z, and colour 





— Infrared - 




P/2MASS £ cdr . 

2MASS* 101 : 4 HiPS for J, H, K, and colour 

9 

0781 

100 % 

P/2MASS6X £ cdr . 

2MASS6X: 3 HiPS for J, H, and K 

9 

0781 

1.32% 

p/ultravista £cds ... 

VISTA Ultra Deep Survey* 111 : 6 HiPS for Y, J, 

12 

071 

0.005% 


H, Ks, NB118, and colour 




P/WISE £cd r . 

WISE* 121 : 5 HiPS for 3.4, 4.6, 12, 22 pm, and colour 

5 

1279 

100 % 

P/WISE/WSSA £ c d . 

WSSA* 131 Diffuse dust emission at 12 pm 

7 

3722 

100 % 

P/ALLWISE £f£/ ''. 

ALLWISE* 141 : 5 HiPS for 3.4, 4.6, 12, 22 pm, and colour 

8 

1761 

100 % 

P/DIRBE Guv . 

DIRBE* 151 : 20 HiPS for bands 1-10 and ZSMA bands 1-10 

3 

5175 

100 % 

p/iris g uvr . 

IRAS-IRIS* 161 : 4 HiPS for 12, 25, 60, and 100 pm 

3 

5175 

100 % 

p/iris g cdr . 

IRAS-IRIS* 161 : HiPS colour composition 

3 

5175 

100 % 

P/GLIMPSE360 £rwxz .. 

GLIMPSE 360* 17,18) : colour 

9 

0781 

3.1% 

P/SPITZER/IRAC " dyx .. 

Spitzer IRAC: 5 HiPS for the 3.6, 4.5, 5.8, 

9 

0781 

1.4% 


8 /tm, and colour bands 




P/SPITZER/MIPS l c " ;rfvi 

Spitzer MIPS l* 221 24 pm band 

8 

1761 

1 .6% 

P/SPITZER/MIPS2 g cdyx 

Spitzer MIPS2* 221 70 pm band 

7 

3722 

1 .6% 

P/SPITZER/MIPS3 g cdyx 

Spitzer MIPS3* 221 160 pm band 

6 

6744 

1 .6% 

p/akari/fis £ cda . 

AKARI* 231 Far infrared all-sky survey: 5 HiPS for N60, 

5 

1279 

99.9% 


N160, WideS, WideL, and colour 
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— Sub-mm and radio - 




P/SCUBA/450em £ c d P .. 

SCUBA 450 (24) gm 

8 

1761 

0.3% 

P/SCUBA/850em £ c d/s .. 

SCUBA 850 (24) pm 

7 

3722 

0.9% 

P/SCUBA/850emi £ cd P .. 

SCUBA 850 (24) pm extended data set 

7 

3722 

0.9% 

P/SCUBA2 £<:,r . 

SCUBA-2 (25) : 2 HiPS for 450 and 850 pm 

9 

0781 

1.8% 

P/BOLOCAM 0 “ v . 

Bolocam Galactic plane survey* 261 (BGPS) 

5 

1279 

2% 

P/WMAP C cds . 

WMAP Year 9 results (27) : 5 HiPS for W, V, Q, Ka, and K 

3 

6(87 

100% 


(n,ii e =64) 




P/PLANCK/HFI c c d e ... 

PLANCK (28) (2013 data release): 6 HiPS for HFI in HFI100, 

3 

1(72 

100% 


143, 217, 353, 545, 857 GHz, and colour (n tUe = 256) 




P/PLANCK/LFI 0 (:df .... 

PLANCK (28) (2013 data release): 4 HiPS for LFI in 1030, 

3 

3(44 

100% 


1044, 1070 GHz, and colour (n tae = 256) 




P/PLANCK/CMB c cde .. 

PLANCK CMB fluctuations map (SMICA inpainted) 

3 

1(72 

100% 


(ft tile —256) 




P/NVSS £cc,f . 

The NRAO VLA Sky Survey (29) 

5 

12.88 

84.8% 

P/CHIPASS c " v . 

CHIPASS (30) 1.4 GHz Radio continuum survey 

3 

5175 

72% 

P/SUMSS £crf . 

Sydney University Molonglo sky survey 131,32 - ) 843MHz 

6 

6744 

25% 

P/DWINGELOO c uv .... 

Dwingeloo 820 MHz continuum survey 33 ' 

3 

5175 

57% 

P/HASLAM408 C 71 ed .... 

Haslam 408 MHz survey: 2 HiPS for original 34 ' 

3 

5175 

100% 


and reprocessed' 35 ' data 




P/VLSSR Ecd < . 

VLA low-frequency sky survey redux 36 ' 

5 

1279 

75% 


(74 MHz continuum) 




P/WENSS £c ' i * r . 

Westerbork northern sky survey' 37 ' (325 MHz) 

5 

12.88 

26.9% 


-Ha- 




P/Finkbeiner G c d A . 

Finkbeiner Ha composite survey' 38 ' 

3 

5175 

100% 

P/SHS £crf/i . 

SuperCOSMOS digitized photographic Her Survey' 39 ' 

10 

074 

9.8% 

P/SHASSA C cd ^ . 

SHASSA southern Her sky survey atlas' 40 ': 4 HiPS for 

4 

25777 

67% 


Her, continuum, continuum-corrected. Her 





and smoothed Her images 




P/VTSS Ccdjr . 

Virginia Tech spectral line survey' 41 ': 2 HiPS for 

3 

5175 

20% 


Her and continuum corrected Ha 





-HI- 




P/GASS 0 cdp . 

Parkes Galactic all sky survey' 42 ' (GASS) («,,/,, =64) 

3 

6(87 

53% 

P/GASS/NH 0 " "P . 

Parkes Galactic all sky survey' 42 ' (GASS) (n,u e = 128) 

3 

3(44 

53% 

P/CGPS/VGPS c " vcrT ... 

Combination CGPS' 43 ' and VGPS' 44 ': 2 HiPS for 

4 

2578 

3.9% 


1420 MHz continuum and HI 





-CO- 




P/CO Gds . 

CO composite survey' 45 ' 

3 

5175 

100% 

P/COHRS £c ' /<5 . 

CO high resolution survey' 46 ' (HARP/JCMT) 

7 

3722 

0.15% 
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Notes. <£) Equatorial (ICRS) coordinate system HiPS ; (G) Galactic coordinate system HiPS ; (a) CDS identifier for the HiPS data set ; (i) Original 
data distributed by SkyView/HEASARC ; (c) HEALPixed by CDS ; <d) HiPS distributed by CDS ; (e) HEALPixed by Laurent Michel (Observatoire 
de Strasbourg) ; m HiPS distributed by XMM-Newton science survey consortium, EC FP7 ARCHES project ; Original data distributed by 
INTEGRAL Science Data Centre, Geneva ; (/,) Original data distributed by Space Research Institute (IKI), Moscow ; (,) Original data distributed 
by Max Planck Institut fur extraterrestrische Physik ; ® Original data distributed by MAST at STScI ; ® Original data distributed by the SLOAN 
Digital Sky Survey ; Original data (c) Axel Mellinger ; (n) Original data (c) CFH - powered by Terapix ; (p) Processed data from CADC ; 
(?) HEALPixed by Daniel Durand (CADC) ; (r) Original data distributed by IPAC/NASA Infrared Science Archive (IRSA) ; (s > Original data 
distributed by the Ultra VISTA consortium ; w Processed data obtained from (12) ; >M> HiPS distributed by Analysis Centre for Extended Data 
(CADE), IRAP Toulouse ; (v) HEALPixed by D. Paradis, CADE, IRAP Toulouse ; M HEALPixed by T. Robitaille ; w Original data from 
Spitzer mission, JPL/NASA ; w Composite from Spitzer legacy programs: GLIMPSE (17,18), SAGE (19), SAGE-SMC (20) and SINGS (21) 
; fe) HiPS distributed by Spitzer mission, JPL/NASA ; (a) Original data distributed by ISAS, JAXA, Japan ; ® Processed data obtained from 
CADC ; (y> HiPS built from JCMT data via CADC archive ; (tf) Original data from LAMBDA ; (E) Original data from PLANCK mission at ESA 
;® Original data from NRAO ; 1,1 Data distributed and maintained by LAMBDA, who cite Jodrell Bank Centre for Astrophyscis as the original 
data source ; <ei Reprocessed data by described in (35) including details of various versions of the 408 MHz Haslam map ; {rt Original data from 
WENSS at ASTRON in the Netherlands ; w Composed by VTSS (Virginia Tech spectral line survey), SHASSA (Southern Ho- sky survey atlas), 
WHAM (Wisconsin Her mapper) ; w Original data courtesy of WFAU, AAO/UKST, PPARC/STFC ; ® Original data from Swarthmore College 
Incorporated and available from NASA/Skyview ; M Original data from Virginia Tech Physics ; ® Original data from CSIRO/ATNF ; (cr) Origi¬ 
nal VGPS data from University of Calgary and CADC ; (t) Original CGPS data from Canadian Galactic plane consortium and available from CADC 


References. (1) Atwood et al. (2009); (2) Cillis & Hartman (20051); (3) Winkler et al. (2003); (4))Krivonos et al. 

(2012); (5) Voges (1992); (6) 

Martin et al. (2005); (7) Ahn et al.j<2012); (8) see CFHTLS in the acknowledgements; (9) Giavalisco et al. (2004 

i; (10) Skrutskie et al. (2006); 

(11) McCracken et al. (2012|; (12) Wright et al. (20101; (13) Meisner & Finkbeiner (2014 

>; (14) Cutri et al. (20131); (15) V 

iville-Deschenes & 

Lagache (2005); (16) Boggess et al. (1992); (17) Benjamin et al. (2003); (18) Churchwell et al. 

(2009); (19) Meixner (2008 

; (20) Gordon et al. 

(20111; (21) Kennicutt et al. (20031; (22) Werner et al. (2004); (23) Murakami et al. (2007); (24) 

Holland et al. (1999); (25) H 

olland et al. (2013); 


(26)]Aguirre et al. 

(2011); (27) 

Bennett et al. ^2013); (28) Planck-Collaboration et al. (2014); (29) Condon et al. (1998); (30) Calabretta et al. 

(12014); (31)|Mauc 

het al. (2003 

); (32)|Murphy et al. (2007); (33) Berkhuijsen (1972); (34) Haslam et al. (1982); (35) Remazeilles et al. (2014) 

(36) Lane et al. (2012); (37) Rengelink et al. (1997); (38) Finkbeiner 

(2003 ); (39) Parker et ak|(2005); (40)|Gaustad et al. 

2001); (41)|Dennison 

et al. (1999); (42) McClure-Griffiths et al. (2009); (43) ((Taylor et al. 

2003); (44) Stil et al. (2006); (45) Dame et al. (2001 

i; (46) Dempsey et al. 


(20l5i: (47)|Di Francescoet al. (2008) 
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Table 4. HiPS catalogues 


ID" 

Data set name and description 

ktile,max 

Sort-key 6 

N(sources) 

B/SIMBAD. 

CDS SIMBAD database* llc 

14 

Literature citations 

~8xl0 6 

I/259/TYC2 

Tycho-2 catalogue* 2,3,41 

6 

BT.VT 

2,539,913 

1/312 . 

PPMX catalogue of positions and proper motions* 5,61 

9 

Bl,B2,Rl,Rs,I 

18,088,919 

1/324 . 

Initial Gaia source list* 7 ' 81 

11 

BJ,RF,G,Grvs 

1,222,598,530 

11/246 . 

2MASS* 9,101 

11 

J,H,K 

470,992,970 

11/293 . 

GLIMPSE* 11,12,131 

10 

f3.6,4.5,5.8,8.0 pm 

104,240,613 

11/294 . 

SDSS-DR7* 141 

9 

u,g,r,i,z 

65,714,108 

11/297 . 

AKARI IRC* 15 ’ 161 

7 

S09 

870,973 

11/298 . 

AKARI FIS* i7) 

6 

S65 

427,071 

11/311/WISE 

WISE* 18,191 

8 

W1 band 

563,921,584 

VI/137. 

Gaia universe model snapshot* 20 ’ 211 (GUMS) 

11 

Distance 

2,143,475,885 


Notes. (a> CDS identifier for the catalogue ;* w The sort-key used to distribute catalogue sources over the different orders of the HiPS. Where the 
sort-key is based on source brightness, we indicate the list of the bandpasses that were combined to define the sort-key ; * cl SIMBAD is a living 
database that is continually updated with new astronomical objects, hence N(sources) is growing. 


References. (1) 

Wenger et al. ()2000); (2) H0g et al. 

2000a); (3) H0g et al. 

(2000b); (4) H0g et al. 

2000c 

); (5) Roser et al. 

(2008 

; (6) Roeser| 

et al. (20081; (7) 

Smart & Nicastro (2014); (8) Smart 

(2013); (9) 

Cutri et al. ( 

2003apf (10) Cutri et al. 

2003b 

); (11) Churchwe 

1 et al. 

(2009); (12) 

Benjamin et al. 

2003 

; (13) 
■; (18) 

Spitzer Science Center] i 

2009); (14) 

Abazajian et al. (2009); (15) Ishihara et al. 

(2010apj(l 6) Ishihara et al. (2010b); 

(17) Polio et al. 

2010 

Polio et al. (2010); (19) 

Wright et a 

. (2010); (20) Robin et al. (2012a); (21) Ro 

tin et al. (2012b) 


Table 5. HiPS three-dimensional data cubes 


ID" 

Data set name and description 

ktile,max @pix,min N(pl(lTl€S ) 

Sky 

(%) 

C/EGRET/DIF. 

EGRET Diffuse Maps* 31 

3 5175 10 

100 % 


30-50, 50-70, 70-100, 100-150, 
150-300, 300-500, 500-1000 MeV, 
1-2 GeV, 2-4 GeV, and 4-10 GeV 


C/HST/GOODS 6 c .... 

GOODS* 41 project HST data (b, v, i,z) 

14 

25.1 mas 

4 


C/HST 6 ". 

HST: NICMOS, WFPC2, and ACS in 

FI 10W, F160W, F225W, F300W, 

F450W, F606W, F625W, F702W, 

F814W, and F850LP 

14 

25.1 mas 

13 

0.044% 

C/CALIFA/V500/DR1 

CALIFA DR1 V500 cubes* 61 

9 

0781 

1877 

0.00009% 

C/C ALIFA/V500/DR2 

CALIFA DR2 V500 cubes* 71 

9 

0781 

1877 

0.00018% 

C/ALLWISE . 

ALLWISE* 21 3.4, 4.6, 12, and 22 pm 

8 

1761 

4 

100 % 

C/IRIS . 

IRAS-IRIS* 51 12, 25, 60, and 100 pm 

3 

5175 

4 

100 % 

C/CGPS . 

Canadian Galactic plane survey* 11 

5 

1279 

272 

3.9% 

C/WMAP/CUBE/9YR 

WMAP W, V, Q, Ka, and K (n tUe =64) 

3 

6(87 

5 

100 % 


Notes. , " > CDS identifier for the HiPS cube ; (hl Processed data from CADC ; (ri HEALPixed by Daniel Durand (CADC) 

References. (l)ITaylor et al. 1(12003 b ; (2)|Cutri et al.|(|2013b; (3)|Cillis & Hartman|12005k; (4)|Giavalisco et al.|d2004b; (5)|Boggess et al.|dl992b; (6) 
|Husemann et al.|(|2013|>; (7)|Garcfa-Benito et al.|i|20l4]i 
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